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Abstract. The main goal of this paper is to present a new algorithm bounding the regularity and 
"alpha" (the lowest degree of existing hypersurface) of a linear system of hypersurfaces (in P n ) passing 
through multiple points in general position. To do the above we formulate and prove new theorem, which 
allows to show non-specialty of linear system by splitting it into non-special (and simpler) systems. As 
qq ■ a result we give new bounds for multiple point Seshadri constants on P 2 . 
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■ 1. Introduction 

In what follows we assume that the ground field K is of characteristic zero. Let n > 2, let d S Z, let 
mi, . . . , m r £ N. By C n (d; mi, . . . , m r ) we denote the linear system of hypersurfaces (in P™ :— P™(K)) of 
degree d passing through r points p\,...,p r in general position with multiplicities at least mi, . . . ,m r . 
The dimension of such system is denoted by diraC n (d; mi, . . . , m r ). Define the virtual dimension of 
L = £ n (d;mi, . . . ,m r ) 
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and the expected dimension of L 

edimL := maxjvdimL, —1}. 

Observe that dimL > edimL. If this inequality is strict then L is called special, non-special otherwise. 
The system L is called non-empty if dimL > 0, empty otherwise. 

The problem of classifying special systems has been widely studied by many authors, especially for 
n = 2 (most of papers included in reference section deal with the case n = 2; there are many others). 
For low multiplicities (i.e. bounded by some constant) the classification is known, it begun with [Hi r 85] , 
where the case mi = • • • = m r < 3 was solved; special systems with (arbitrary) multiplicities up to 
11 has been classified in [Dum-Jar 07] ; the homogeneous case (mi = ••• = m r ) with multiplicities 
up to 42 in [Dum 07] : quasi-homogeneous case (mi = ••• = m r _i, m r arbitrary) for mi = 5 (for 
qq ' mi < 5 the classification is also done) in |Laf-Uga 03| . For r < 9 the classical result is due to Nagata 

( Nag 59 ). In the general case (n > 3) we note that very few is known, results concern either low 
multiplicities (mi = • • • = m r = 2, n arbitrary) (see |Ale-Hir 95j . [Cha 01| ). or few points (r < 8, n = 3 
in |Vol-Laf 07] V 

For a system of multiplicities mi, . . . , m r and n > 2 we define (after Harbourne, [Har 02j ) 

H [ 

a n {mi, . . . ,m r ) := min{d G N : £ n (d;mi, . . . ,m r ) is non-empty}, 

T n (mi, ■ ■ ■ , m r ) :— minjd € N : C n {d; mi, . . . , m r ) is non-special and non-empty}. 

The second number is often called the regularity. Observe that for d < a n (mi, . . . ,m r ) the system 
C n {d;mi, . . . ,m r ) is empty, similarly, for d > r„(mi, . . . , m r ) the system C n (d; mi, . . . , m r ) is non- 
special. More information about ct2 and Ti can be found in |Har 02] and [Har 05j . Note that usually (for 
small values of n) the conjectured values of a n and T n are known (for n = 2 by Harbourne-Hirschowitz 
conjecture (see e.g. [Hir 89] ): for n = 3 the conjecture can be found in |Laf-Uga 06| ; for n > 4 in many 
cases we can guess), which gives upper bound for a n and lower bound for t„. In this paper we present a 
new algorithm to bound a n (from below) and t„ (from above). We study the behavior of this algorithm 
for T2 and a n for n — 2,3,4. Additionally we give new bounds on multiple point Seshadri constants on 
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The paper is organized as follows: Next section is devoted to present main facts used in the paper. 
Theorem [T] is the main tool and (as author believes) is interesting on its own. Other theorems are (more 
or less) known. In section 3 we present two algorithms (called NSsplit and NSglue). Both are used 
to show non-specialty of a given system. The first one tries to split system into many smaller (i.e. lower 
degree, less multiplicities) systems (it is more probable to find a suitable criterion of non-specialty for 
a small system; some of such criterions are given in section 2). The second one uses Thm. Q]to „glue" 
multiplicities, which allows usage of birational isomorphism (Thm. [3]). In section 4 we present examples 
of bounds on T2, a^, 0:3 and 0:4. We focus on quasi- homogeneous systems (for explanation see Thm. 
H"4"|) . The bounds (especially for regularity) are (as far as author knows) much better than bounds known 
before. In the last section we show how our method, together with Thm. [16] (proposed by Eckl), can be 
used to produce bounds on multiple point Seshadri constant on P 2 . 



2. Main tools 

In |Cil-Mir"98] Ciliberto and Miranda proposed a method of computing the dimension of a linear 
system in P 2 by splitting the problem into several ones (possibly easier). Similar method (based on 
splitting) was used by Biran ( |Bir 99] ) to prove ampleness and nefness of divisors on blow-ups of P 2 . 
Our theorem uses the same concept, but the proof is completely different and works in any dimension 
(n > 2). 

Theorem 1. Let n > 2, let d, k, mi, . . . , m r , m r +i, . . . , m s S N. If 

• L± = C n (k; mi, . . . , m s ) is non-special, 

• L% = C n {d; m s _|_i, . . . , m r , k + 1) is non-special, 

• (vdimLi + l)(vdimL 2 + 1) > 0, 

then the system L = C n (d; mi, . . . , m r ) is non-special. 

In the proof of Theorem [T] we will use the reduction method introduced in |Dum 07] , We must adapt 
this method to the n-dimensional case. The following notations are used only in the proof of Theorem [T] 

Every system C n (d; mi, . . . , m r ) can be considered as a vector space V n (d; mi, . . . , m r ) (over K) of 
polynomials / G K[x] = K[xi, . . . , x n ] of degree at most d, such that mult Pj . / > rrij, for p\, . . . ,p r € K™ 
in general position. We have 

diniK V n (d; mi, . . . , m r ) = dim£ n (d; mi, . . . , m r ) + 1. 

Let 

vdim V n (d; mi, . . . , m r ) := vdim C n (d; mi, . . . , m r ) + 1, 
cdim V n (d; mi, . . . , m r ) := cdim£„(d; mi, . . . , m r ) + 1. 

The space V — V n (d; mi, . . . , m r ) is said to be non-special if dimu V = cdim^, which is equivalent to 
the non-specialty of C n (d; mi, . . . , m r ). 

By V n (d) we understand the space of all polynomials of degree at most d with no conditions imposed. 
Let 

Vj : V n (d) ^kC" 3 """ 1 ) 

be a linear function which maps / to a collection of all partial derivatives of / up to degree rrij — 1 
evaluated at pj . Put 

<S>:V n (d)3 f^( Vl (f ),..., Vr (f))£K c , 

where 

^m., + n — 1 N 
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It is easy to observe that V n (d; mi , . . . , m r ) — ker <E>. Let M n (d; mi , . . . , m r ) be the matrix of $ in mono- 
mial basis of V n (d) and canonical basis of K c . Columns of M n (d; mi, . . . , m r ) are indexed by monomials 
of degree at most d, rows are indexed by conditions (i.e. points and symbols of partial derivatives). This 
matrix depends on coordinates of points pi, . . . ,p r , hence we consider each entry of M n (d; mi, . . . , m r ) 
as a polynomial (in fact, it is a monomial with coefficient) in P := K[p}, . . . ,p" , . . . ,pj,, . . . where 
Pj is the fc-th coordinate of j-th point. Observe that system C n (d; mi, . . . , m r ) is non-special if and only 
if the matrix M n {d;m\, . . . ,m r ) has maximal rank (as a matrix over P; in other words, there exists a 
square submatrix M of M n (d; mi, . . . , m r ) of maximal size such that det M is a non-zero polynomial in 

py 
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Define the multidegree function mdeg on P by setting 

mdeg(^) := (0, . . . , 0, 1, 0, . . . , 0) G N™ 

with 1 on fc-th position. 

Lemma 2. If M is a square submatrix of M n (d; mi, ■ ■ ■ , m r ) then detM is multi-homogeneous with 
respect to mdeg. Moreover, 

mdeg(det M) = - 7, 

where (3 is equal to the standard multidegree (in K[xi, . . . , x n \) of product of monomials indexing columns 
of M, and 7 depends only on conditions indexing rows of M . 

Proof. The proof can be done easily by examining the entries of M, and is esentially the same as the 
proof of Proposition 11 in [Dum 07] , □ 

Proof of Theorem\]$ Consider the space V n (d; m s +i, ■ ■ ■ , m r , k + 1). Since a translation T : K™ — ► 
K™ produces an isomorphism of V n (d; m 3 +i, ■ ■ ■ , m r , k + 1) based on p s +\, ■ ■ ■ ,Pr,Pr+i with the same 
space based on T(p s+ i), . . . ,T(p r ),T(p r+ i), we may assume that p r +i = (0, ...,0). Now the space 
V n (d; m s +i, . . . , m r , fe + 1) is nothing else than the space V n (d, k; m s +i, . . . , to,.) of polynomials generated 
(over K) by the set 

{x 13 : k < deg x < d} 

with respective multiplicities at p s +\, ■ ■ ■ ,p r - 

Let us consider the matrix M = M n (d; mi, . . . , m r ). If necessary, we may reorder columns and rows of 
M such that the first ( fe ^ n ) columns are indexed by monomials of degree at most k, first 53^=1 ( mj ^ l ™ _1 ) 
rows indexed by points pi, . . . ,p s . Let Mi be the left upper submatrix of M consisting of ( fe ^ n ) columns 
and EU C^tr 1 ) ^ 



M = 



Mi 



M 2 



K 2 

Observe that Mi = M n (k; mi, . . . , m s ), M 2 — M n (d, k; m s +i, . . . , m r ), so both Mi and M% have maximal 
rank. 

Now, let us consider the case vdimLi > —1 and vdimL2 > — 1. Then rows(Afi) = rank(Mi) < 
cols(Mi), rows(M 2 ) = rank(M 2 ) < cols(M 2 ). Take a square submatrix M[ of Mi (resp. M 2 of M 2 ) such 
that rows(Af 1 ) = rows(Mi), detM{ ^ (resp. rows(M2) = rows(M 2 ), detM 2 7^ 0). Consider the matrix 

M' = 



M[ 




K' 2 


M' 2 



Observe that M' is a square submatrix of M of size rows(M). It is enough to show that det M' 7^ to 
complete the proof. 

Let Di be the set of monomials indexing columns of ML I — 1,2, put D := DiUD 2 . Let U = [M[ \ K[], 
L = [K' 2 I M 2 ] be submatrices of M' , let 

C:={CcD|#C=Pi}. 

For C C D define Uc (resp. Lc) as the submatrix of U (resp. L) given by taking the columns indexed 
by elements of C. Now we can compute 

det M e ( C ) det U C det L D \ C , 

cec 

for e(C) = ±1. Observe that 

detUceK\ P 1 1 ,...,pl,...,pl,...,p*}, 
det L D \ C G K[pl +1 , . . . ,p™ +1 , . . . ,pl, . . . ,p™\. 

We have 

det M' = ± det M[ det M 2 + f. 
Assume that det M' = 0. As detMj ^ and detM 2 ^ 0, there exists C £ C, C / Di, such that 
mdeg(detC/ c ) = mdeg(detM{). By Lemma [2] 

mdeg(det U c ) = ( ^ 0) - 7, mdeg(det M[) - ( ]T 0) - % 
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so 

E ifl = E ifl- 

However, for any a;^ 1 S Di and a;' 32 £ D \ D\ we have < j/^l, which leads to a contradiction. The 
case vdimLi < — 1 and vdimL2 < — 1 can be done similarly. □ 

Theorem 3. Let n > 2, let d, mi, . . . , m r € N, let k = (n — \)d — (mi + • ■ ■ + m„+i). If rnj + k > for 
j = 1, . . . , n + 1 i/ien 

dim£„(<i; mi, . . . , m r ) = dim Cr(£ n (d; mi, . . . , m r )), 

where 

Cr(£ n (d;mi, . . . ,m r )) := £„(d + k;mi + k, ... ,m n +i + k,m n+2 , ■ ■ ■ ,m r ). 

Proof. The proof for n = 2 is well-known. In Laf-Uga 06 we can find the proof for n — 3 and that 
idea can be applied to arbitrary n. Namely, using projective change of coordinates assume pi, ■ ■ ■ ,p n +i 
to be fundamental ones. The system £ n (d;mi, . . . , m r ) is equivalent to the system of hypersurfaces 
C„ (M; m n+ 2, ■ ■ ■ , m r ) generated by the set 

M = {xf 1 • . . . • :\[3\=d,{3 j <d-m j ,j = l,...,n+l}. 

Similarly, system L n (d + k; mi + k, . . . , m n +i + k, m n +2, . . . , m r ) is equivalent to £ n (M'; m„ + 2, . . . , m r ) 
for 

M' = {x^ 1 x^ 1 :\j3\=d + k, p 3 <(d+k)- (rrij + *), j = 1, . . . , n + 1}. 
The standard birational transformation 

P* : {xi : . . . : x n+1 ) .— > (.t^ 1 : . . . : z^) € P n 

induces a bijection 

M 3 X? 1 ■ . . . ■ X^ .— > X d - mi - Pl ■ . . . • x d n -^ + l-0n+l g M / ; 

which completes the proof. □ 

Geometrically speaking, the system Cr(L) is an image of L by a birational transformation. Such 
transformation (in P 2 ) is often reffered as to Cremona transformation. In what follows we use the name 
"birational transformation" to denote Cr operation. 

Theorem 4. Let n > 2, let d, mi, . . . , m r G Z. If (n—l)d— 2j=i m j < nij > /or j = 1, . . . , n i/ien 

dim C n {d\ mi, . . . , m r ) = dimHp(£ n (<i; mi, . . . , m r )), 

w/iere 

Hp(£ n (d;mi, . . . ,m r )) := C n (d - l;mi - l,...,m„ - l,m n+ i, . . . ,m r ). 

Proof. The proof for n = 3 can be found in |Laf-Uga 06] , we use the same reasoning. Consider the 
situation as in the previous proof. Since Pi + ■ ■ ■ + f3 n < nd — (mi + • • • + m„) < d, the map 

IVl d>\ ... X n+ l < * X\ ■ ■ ■ a„ ■i„_|_l 

is one-to-one and its image corresponds to the system C n {d — 1; mi — 1, . . . , m n — 1, m n +i, . . . , m r ). 
Geometrically speaking, the hyperplane passing through the first n points is in the base locus of 
C n (d;mi, . . . ,m,) and can be "taken out". □ 

Definition 5. We say that C n (d; mi, . . . , m r ) is in standard form if d < or the following holds: 

• mi , . . . , m r are non- increasing, 

• (n-lJi-Egnij > 0. 

Proposition 6. i?ac/i system L has its standard form L' such that dimL = dimL'. 

Proof. The following set of instructions allows to transform any system L into a standard form. 

• sort multiplicities in non-increasing order, 

• if (n — l)d — Y^j=x m j < then take L < — Hp(L) and go back to the first step, 

• if (n — l)d — J2j=i m j < then take L < — Cr(L) and go back to the first step. 

Observe that each time (in steps 2 and 3) the degree of L decreases and the dimension does not change. □ 
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The end of this section is devoted to prepare a family of non-special systems (which will be used later) 
by presenting criterions of non-specialty. We will use the following notation: m xk denotes the sequence 
of m's taken k times, 

m xk = (to, . . . , to) . 

k 

Criterion 1. Let C n (d; m\, . . . , m r ) be a standard form of a system L. If either d < or mi < 1 then 
L is non-special. 

Proof. The case d < is obvious. If mi < 1 and C n (d; mi, . . . , m r ) is in standard form then rrij < 1 for 
j = 2, .. . ,r. Each point of muliplicity one (or zero) imposes an independent condition. □ 

Remark 7. It is known that we can perform Cr operation on the system with (n — l)d — Y)j—i mj < 0. 
This leads to negative multiplicities, which can be properly understood using divisors on the blow-up of 
P™ in r points. However, it is easy to show (using Thms. [3] and [3]) that if L = C n (d; mi, . . . , m r ) can be 
transformed (by birational transformations) to C n (d'; to' 1; . . . , m! r ), for m\, . . . , m' r € Z and d < 0, then 
L is empty. 

Criterion 2. Let L — Li(d;m\, . . . ,m r ). If vdimL > —1 and L can be reduced to a system L' (using 
reductions decribed in [Dum-Jar 07] ; the system L' can be generated by some finite set of monomials) 
such that vdimL' = vdimL and L' has no conditions (multiplicities imposed) then L is non-special. 

Proof. By Thm. 17 in [Dum-Jar 07| dimL < dimL'. Since L' is condition-free, we have dimi' = 
vdimL' = vdimL. □ 

Criterion 3. Let L = £2^; "ii, ■ ■ ■ ,m r ). If vdimL < —1 and L can be reduced to a system L' (using 
weak reductions decribed in [Dum-Jar 07] ) such that L 1 has no monomials, then L is empty (and hence 
non-special). 

Proof. Again, by Thm. 17 in [Dum-Jar 07| we have dimL < dimL' = — 1. □ 
Criterion 4. Let to > 0. The systems Lii^m — 1; to x4 ) and Lii^m; to x4 ) are non-special. 
Proof. Use birational transformation. 

Criterion 5. Let r > 3, let to > (r - 2)/4. Then r 2 (m xr2 ) = rm + \(r - 3)/2] 
Proof. See [HHF 03) . Lemma 5.3. 

3. Algorithms 

To bound a n and r„ it is sufficient to find an algorithm, which, given n, d, m±, . . . , m r , returns either 
non-special (and then the system C n {d;m\, . . . ,m r ) is non-special), or not-decided. We will focus 
on such algorithms. The answer SPECIAL is also allowed (of course only if C n (d; m\, . . . , m r ) is special). 

Let us assume we are given a collection (called Eclass) of non-special systems. We will assume that 
if L e Eclass then also Cr(L) and Hp(i) are in Eclass. The first algorithm (NSsplit) makes use of 
Theorem [TJ and tries to show non-specialty of a given system by splitting it into systems belonging to 
Eclass. 

Algorithm NSsplit 

Input: n > 2, d, mi, . . . , m r € N 
Output: NON-SPECIAL or NOT-DECIDED or SPECIAL. 

vi < — vdim C n (d; toi , . . . , m r ) ; 
change C n {d; m\, . . . , m r ) into standard form; 
if vdim C n {d; mi, . . . , m r ) > max{0, v\] then return SPECIAL; 
if £ n (d;mi, . . . , m r ) £ Eclass then return non-special; 
for each \i\, . . . , i s } C {1, . . . ,r} and < k < d do 
{h,---,3r-s} < — {l,...,r}\{£i,...,i 8 }; 

v < — (vdim£ n (/c; to^ , . . . , to^ s ) + l)(vdim C n (d;k + 1, , ■ ■ • , mj r _„) + 1); 

Al < — NSsPLlT(n, k, TOjj, . . . , TO Js ); 

A2 < — NSSPLlT(n, d, k + l,m jl , m. ]r _ s ); 

if v > and Al = non-special and A2 = non-special then return non-special; 



□ 
□ 



end for each 
return not-decided; 

The above algorithm can be implemented in a more subtle way, for example for v < we can skip the 
next two steps (running NSsplit). Also one can use criterions of specialty to avoid running NSsplit 
on special systems. 

Example 8. We will show that £2(72; 10 x50 ) is non-special (i.e. empty). Applying NSsplit with 
ECLASS := {L = £ 2 (d;m 1 , . . . , m r ) : L is empty by Crit. Q]} 

we can find the following solution (presented in 5 steps): 

(1) split £ 2 (72;10 x50 ) into £ 2 (72; 61, 10 x15 ), which is empty by Crit. [fl and £ 2 (60; 10 x35 ); 

(2) split £ 2 (60; 10 x35 ) into £ 2 (45; 10 x20 ) and £ 2 (60;46, 10 x15 ), which, by Cremona transformation, 
is equivalent to the system £ 2 (18; 10,4 x15 ); 

(3) split £ 2 (18; 10,4 x15 ) into £ 2 (10;4 x7 ) and £ 2 (18; 11, 10, 4 xS ); both systems are empty by Crit. 

m 

(4) split £ 2 (45; 10 x20 ) into £ 2 (28; 10 x8 ) and £ 2 (45; 29, 10 x12 ); the first system is empty by Crit. Q] 
the second one can be transformed into £ 2 (21; 6 xl5 , 5); 

(5) split £ 2 (21;6 xl5 ,5) into £ 2 (11; 6 x4 ) and £ 2 (21; 12, 6 x8 , 5); both are empty by Crit. HJ 

Algorithm NSsplit is very slow for systems with many multiplicities. Therefore we will present 
another algorithm, which does not search through the tree of all possibilities, but tries to "glue" some 
multiplicities. The aim is to obtain multiplicities relatively big with respect to the degree, i.e. (n— l)d < 
Y^j—y Trij, and use birational transformation(s). 

Theorem 9. Let £„(fc;m s ) be non-special, let 

Li = C n (d; mi, . . . , m r , m s ), 

L 2 = £ n (d; mi, . . . , m r , k + 1). 

// either — 1 < vdimL 2 < vdimLi or vdimLi < vdimL 2 < — 1 then in order to show non- specialty of 
L\ it is enough to show non-specialty of Li. 

Proof. Follows from Theorem [TJ □ 

The algorithm NSglue depends on two sets (ECLASS and Gclass) of non-special systems. The set 
GCLASS containing systems of the form C n (k; m xs ) should be relatively small and should contain systems 
with virtual dimensions close to —1. 

Algorithm NSglue 

Input: n > 2, d, mi, . . . , m r e N 
Output: NON-SPECIAL or not-decided or SPECIAL. 

v\ < — vdim £„ (d; mi , . . . , m r ) ; 

change £ n (d; mi, . . . , m r ) into standard form; 

if vdim£ n (<i; mi, . . . , m r ) > max{0, vi} then return SPECIAL; 

if £ n (d; n%i, . . . , m r ) e Eclass then return non-special; 

V2 < — vdim£„((i; mi, . . . , m r ); 

for each £ n (fc; m xs ) 6 GCLASS such that m^ = • • • = m.; s = m for {iy, . . . , i s } C {1, . . . , r} do 

{jl, ■ ■ -Jr-s} < {1, . . . , r} \ {i 1 , i s }; 

V3 < — vdim£ n (d; k + 1, irij 1 , . . . , rrij r _ s ); 
A < — NSGLUE(n, d, k + 1, rrij 1 , . . . , mj r _J; 

if (—1 < v 3 < v 2 or v 2 < v 3 < —1) and A = non-special then return non-special; 
end for each 
return not-decided; 

Observe that vdim£„(d; mi, . . . , m r , fc+1) = vdim£„(d; mi, . . . , m r , m xs ) — vdim£„(fc; m xs ). There- 
fore in order to show non-specialty of a system with non-negative (resp. non-positive) virtual dimension 
we must use "glueing" systems (from Gclass) with non-negative (resp. non-positive) virtual dimension. 
In both cases the virtual dimension after "glueing" will go closer to —1. 





Example 10. Let us show how NSglue works with the system £ 2 (5918; 4000, 1000 x19 ). Put 

ECLASS :— {L : L is non-special by Crit. H]or[2][, 

GCLASS := {£ 2 (2m;m x4 )}. 

By Crit. H]the set GCLASS contains non-special systems of dimension m + 1. 

The standard form of £ 2 (5918; 4000, 1000 x19 ) is £ 2 (5180; 3262, 1000, 918 x18 ). Now we glue the four 
points (918 x4 — > 2 • 918+ 1 = 1837) obtaining £ 2 (5180; 3262, 1837, 1000, 918 x14 ), which can be trans- 
formed into its standard form £ 2 (4261; 2343, 918 x15 , 81). Again, we glue the points to consider the system 
£ 2 (4261; 2343, 1837, 918 x11 , 81). This can be transformed into £ 2 (3424; 1506, 1000, 918 xl °, 81 x2 ). The 
last glueing gives £ 2 (3424; 1837, 1506, 1000, 918 x6 , 81 x2 ). The standard form is £ 2 (112; 22, 9 x3 , 4 x7 ). By 
Crit. [2]the last system is non-special. Alternatively, we can show that r 2 (22, 9 x3 , 4 x7 ) > 31. Glueing 4 
points of multiplicity 4 we get the system £ 2 (31; 22, 9 x4 , 4 x3 ). The standard form is equal to £ 2 (13; 4 x4 ), 
which is non-special. 

Example 11. Let us consider the system L = £ 2 (<i; m x21 , h). We will show that if 197d— 42(21m + h) < 
then L is empty. By Thm. [5]it is enough to show that L' = £ 2 (d; (2m) x4 , m x5 , h) is empty The system 
L' has 10 base points, say pi,. . . ,pi with multiplicity 2m, p^, . . . ,pg with multiplicity m and pio with 
multiplicity h. We perform birational transformations based on triples of points (allowing negative 
multiplicities, see Rem. [TJ, according to the following sequence (of numbers of points): (1, 2, 3), (4, 5, 6), 
(4,7,8), (4,9,10), (1,2,3), (5,6,7), (8,9,10), (5,6,7), (1,2,3), (4,5,6), (4, 7,8), (4,9 ,10), (1,2,3). After 
all computations (done by hand or with Singular procedure available at [Dum 08] ) we get the system 
£ 2 (197d - 882m - 42h; (84d - 376m - 18h) x4 , (42d - 188m - 9h) x 5 , A2d - 189m - 8h). This result allows 
to bound the Seshadri constant of 0^(1) for 22 points in general position (see Prop. fTS"]) . 

4. Results for n = 2, 3, 4 

4.1. Case n = 2. For n — 2 many algorithms bounding a 2 and r 2 are known (see e.g. [Har 03) . |Roe Ola] 
and |Roe Olbj for Roe's unloading method, [Har-Roe 04] for modified unloading, |Mon 07] . and finally 
Har 02] for survey through many algorithms). We will use algorithm NSglue with 

ECLASS = {L : L is non-special by Criterion [TJ [5] or [3]} , 

GCLASS = {£ 2 (2m; m x4 )} U {£ 2 ((2s + l)m + s - 1; m x ^ s+1 '> 2 ), s = 1, . . . , 10}. 

By Crit. [5] the set GCLASS contains non-special systems (this criterion does not work for all m, but 
even in the worst case it works for m > 5; for m < 4 everything is known, see e.g. [Cil-Mir 00] ). The 
implementation of NSglue with the above classes can be found and downloaded from [Dum 08 j. The 
work is done in FreePascal. 

For a very large family of quasi-homogeneous systems our algorithm gives better bounds than any 
other one. In [Har Roc 07J we can find a list of homogeneous systems for r = 10, . . . , 99 being not a 
square (one system for each r). They are conjectured to be empty, but the authors of [Har-Roe 07] were 
unable to show this using known algorithms. For 64 of these systems the algorithm NSglue gives the 
expected answer, for 18 we must use NSsplit. Only for r = 10, 11 our algorithms do not work. 

We will present only few examples, first two of them appearing in [Har-Roe 07] . four of them proposed 
in |Mon 07] , In the Tab. Q] one can find results the auhor was able to find using known algorithms 
(for convenience, we present the difference between bound and conjectured value). For a 2 they were 
all obtained by one of Harbourne/Roe's algorithms (using [Har 08] or by author own implementation 
for quasi-homogeneous case; the author is not sure that they are the best), for t 2 the bounds were 
either computed with |Har 08] (homogeneous case), or obtained by Monserrat (quasi- homogeneous case; 
[Mon 07] ). who checked other algorithms and claimed to find the best bounds. 

Following Harbourne f [Har 02j ) we present the results of NSglue for homogeneous systems with 
1 < m < 150, 10 < r < 450. On the graph presented on Fig. [T] each dot denotes the success — the 
bound of r 2 is equal to its conjectured value, the same rule is used to present results for a 2 on Fig. 
[2] The graphs presenting partial successes (the difference is not greater that some constant) and those 
situations, when bound for a 2 is not less that conjectured by Nagata (see the last section) can be found 
at |Dum 08j . 

Bounding r 2 with NSglue is much easier than bounding a 2 , since, for any m, there exists non-special 
system with very low positive virtual dimension (vdim£ 2 (3m; m x9 ) = 1), which is not true for systems 
based on s 2 points and negative virtual dimension. 
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system 


conj. a 2 = t 2 


Har-Roe a 2 


Mon/Har t 2 


NSglue a 2 


NSglue t 2 


(24 x1 ^) 


84 


-1 


+1 





+1 


(173 xyti ) 


1699 


-19 


+4 


-3 


+1 


(4000,1000 xiy ) 


5917 


-255 


+92 


-1 


+1 


(6000,1500 xiy ) 


8875 


-382 


+138 


-1 


+1 


(500 xiUUU ) 


15826 


-25 


+ 16 


-3 


+0 


(1200 X1UUU ) 


37962 


-41 


+ 17 


-6 


+0 



Table 1 . Examples of bounds for a 2 and r 2 



A 




number of points r € [10, 450] 



Figure 1. Graph showing when NSglue gives conjectured value of r 2 



A 




number of points r <G [10, 450] 



Figure 2. Graph showing when NSglue gives conjectured value of a 2 



4.2. Case n = 3,4. For n — 3 and n — 4 very few is known. The systems with multiplicities all equal 
to 2 have been characterized in [Ale-Hir 95] . For r < 8, n = 3 the necessary and sufficient condition of 
non-specialty can be found in [Vol -Laf 07 . For r < n + 1 the methods presented in the proof of Thm. [3] 
are sufficient to decide about non-specialty. For r > 8 there are, so far, no nontrivial algorithm to bound 
c*3 and «4. 
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For n = 3we will consider homogeneous systems and use NSglue with 
ECLASS = {L : L is non-special by Crit. [J}-, 

GCLASS = {L = C 3 (d; m xr ) : r G {4, . . . , 8} and L is non-special by Crit. [T]}. 

Example 12. We will show that a3(5 x50 ) = 20. Since (by computation of virtual dimension) we have 
a 3 (5 x5 °) < 20, it is enough to show that £3(19; 5 x5 °) is empty. Observe that (by Crit. [T]) systems 
£3(9; 5 x7 ) and £3(8; 5 x5 ) are empty. By "glueing" points (six times using £3(9; 5 x7 ) and then £3(8; 5 x5 ) 
once) it is enough to show that £3(19; 10 x6 , 9, 5 x3 ) is empty. This is true by Crit. [T] 

Some bounds for a3(m xr ) obtained with NSglue are presented in Tab. [2] First we present the 
conjectured value, then (after comma) our bound. 



n = 3 


r = 30 


r = 50 


r = 100 


r = 150 


777 = 5 


17,16 


20,20 


26,24 


30,28 


777 = 10 


33,31 


39,37 


49,46 


57,52 


777 = 20 


64,61 


76,74 


96,91 


110,104 


777 = 50 


157,151 


186, 182 


235,226 


269, 258 


777 = 100 


312,301 


371,364 


467,452 


535,514 



Table 2. Some bounds for a 3 obtained with NSglue 



The same can be done for n = 4 (here, in Tab. we present results bounding 04, then we give an 
example). 



77 = 4 


r = 30 


r = 50 


r = 100 


777 = 5 


13,12 


15,14 


18,16 


777 = 10 


25,23 


28,26 


34,31 


777 = 20 


48,44 


55,51 


66,59 


777 = 50 


119,110 


135, 126 


161,147 



Table 3. Some bounds for a 4 obtained with NSglue 



Example 13. The conjectured value of a 4 (5 x50 ) is 15. We will show that a4(5 x50 ) > 14. By Crit. [T]the 
system £4(8; 5 x8 ) is empty. We begin with system £4(13; 5 x5 °), glue points to obtain £4(13; 9 x6 , 5 x2 ) 
and use Crit. [T]to show that the last system is empty. 

Our algorithm can also be used in higher dimension. Observe that usually, for 77 > 3, if (77 — l)d — 
m j < then vdim Cr(£„(d; mi, m r )) > vdim£ n (d; mi, . . . , m r ). It follows that the birational 
transformation is useless in proving regularity for 77 > 3. 

5. Bounds for Seshadri constants 

Let pi , . . . , p r £ P 2 be points in general position. Define the Seshadri constant of the line bundle 
CV(1) (see jDem 92j or jSze 01] ) 

£(C P 2 (1), r) := inf ( ^ r deg ^ : C is a curve on P 2 

I Ej=i mult w C 

The famous Nagata Conjecture (see |Nag 59| ) states (in the language of Seshadri constants) that 

e(P pa (l),r) = 4= 

for r > 9. The conjecture is still open. It is known that e(Op2(l),r) < -^=, so we may look for lower 
bounds for Seshadri constant. We will use the following fact: 

Theorem 14. Any irreducible and reduced curve C C P 2 disproving Nagata Conjecture for r > 10 is 
almost homogeneous, that is there exists d,m, tuq € N such that C € £2(^5 m xr ~ 1 , mo). 

9 



Proof. See |Syz 05] , Prop. 2.10, where the above theorem is stated for a polarized surface (X, L) with 
Picard number equal one, two or three, or |Har-Roe 07] . Prop. 3.2, where additional conditions on 
d, m, mo are given. □ 

Proposition 15. 

£ (O p2 (l),22)>^ = 



22 



l 

1764 



Proof. Observe that it is enough to exclude the existence of an irreducible and reduced curve C C 
such that 

degC* 42 



T,j=i mult 



Pi 



C < 197' 



42 
197 



By Thm. [Ml such curve C would be a member of £2(6?; m x21 , h) for some d,m, h G N such that 
The last system is empty (see Ex. ITT]) . 

Using the same technique one can show that (see |Dum 08] ) 

39 



21m+h 



< 
□ 



e (O p2 (l),19) > — = 



19 



1 

1521 



For r = 10 the above technique gives the lower bound equal to ^ (glue £2(^5 m x9 , h) to £2(^5 2m, m x5 , h); 
consider cases m < h and m > h separately), but much better bounds are known (see [Har-Roe 07] ). 
We note here that for r = 19 and r = 22 our bounds are equal to bounds presented in [Har-Roe 07j 
(given after Biran |Bir 99] ). 

Since NSglue works better for regularity (than for 0(2) we will rather make use of the following 
Theorem (proposed and proved by Eckl in [Eckl 05| ): 

Theorem 16 (Eckl). Let r > 9 ; let (dj,rrij) be a sequence such that 
• dj >T 2 ((m J + l) x a 

. 4_ 4 > 1. 

Then 

e(Op»(l),r)>a4=- 

To estimate the Seshadri constant for r points we run NSglue to find d(r,m) > T2(m xr ) for m 
varying from 1 to some M € N, depending on time of computations. For each m and d(r, m) we consider 

satisfying the assumptions of Eckl's 

- 1 



d{r,m) 

Theorem. Hence 



yfr and constant sequence (dj,rrij) — (d(r, m),m — 1) 



e(0p>(l),r) > 



d(r, m) 



d(r, m) 



The bound is given by the maximum of 

m — 1 



m = 1, 



,Af. 



d(r, m) 

A suitable computer program can be found at |Dum 08"] . Here (Tab. [¥]) we present our bounds for 
r = 10, . . . , 24 together with non-special system giving the above maximum. Note that for r = 15, 20, 24 
our bounds are slightly better than bounds presented in [Har-Roe 07] , 



r 


bound 


system 


r 


bound 


system 


10 


100/331 


£2 (26480; 8001 xiU ) 


18 


8096/34413 


£ 2 (34413;8097 xiS ) 


11 


8098/26997 


£ 2 (26997;8099 xli ) 


19 


39/170 




12 


8098/28347 


£ 2 (28347;8099 xi ^) 


20 


1617/7235 


£ 2 (36175;8086 x;2U ) 


13 


8098/29697 


£ 2 (29697;8099 xia ) 


21 


4000/18403 


£ 2 (36806;8001 x:ii ) 


14 


4049/15186 


£ 2 (30372;8099 x14 ) 


22 


42/197 




15 


426/1651 


£ 2 (31369;8095 xl! ') 


23 


8095/38862 


£ 2 (38862;8096 x;ia ) 


17 


2000/8473 


£2(33892;8001 X1Y ) 


24 


8092/39657 


£ 2 (39657;8093 XM ) 



Table 4. Lower bounds for e(C P 2(l), r) for r — 10, . . . , 24 
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Studying behavior of NSglue for n — 2 and homogeneous systems we may observe that it works very 
well for r = k 2 — 1, on the other hand the worst results are obtained for r = k 2 + 1 and, surprisigly, 
r = 13. 

One may ask if NSglue can compute Seshadri constant, i.e. if 

in — 1 

lim = e(Opa(l), r), 

m— »oo a(r, m) 

where d{r, m) is the bound on T2(m xr ) obtained with our algorithm. The answer is not known. 
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